Modeling gravitational recoil from precessing highly-spinning unequal-mass black-hole 

binaries 
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We measure the gravitational recoil for unequal-mass-black-hole-binary mergers, with the larger 
BH having spin a/m H — 0.8, and the smaller BH non-spinning. We choose our configurations such 
that, initially, the spins lie on the orbital plane. The spin and orbital plane precess significantly, and 
we find that the out-of plane recoil (i.e. the recoil perpendicular to the orbital plane around merger) 
varies as n 2 /(1 + q), in agreement with our previous prediction, based on the post-Newtonian scaling. 
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I. INTRODUCTION 

The recent observational discovery of a possible recoil- 
ing supermassive black hole at a speed of 2650 kms -1 
with respect to its host galaxy [l[ represent the first ob- 
servational evidence in support of the predictions of Gen- 
eral Relativity in the strong-field, highly-dynamical, and 
highly-nonlinear regime. This recoil, if it in fact resulted 
from a black-hole merger, would confirm the theoretical 
prediction of Campanelli et al. [2J that black-hole mergers 
can lead to very large recoils. This original prediction, 
and the subsequent calculation Q that indicated that 
such recoils can be as large as 4000 kms -1 , has been 
a key trigger for the recent theoretical efforts to deter- 
mine EM signatures of large recoils and the astronomical 
searches for these signatures [E 5], IE IE S Hi • 

Thanks to recent breakthroughs in the full non-linear 
numerical evolution of black-hole-binary spacetimes [9j, 
ITol . [TTj |. it is now possible to accurately simulate the 
merger process and examine its effects in this highly non- 

lincar regime [H EE El EE EE EE EE EE Hi HE Hi m, 

I24I I25I [2q | . Black- hole binaries will radiate between 2% 
and 8% of their total mass and up to 40% of their angu- 
lar momenta, depending on the magnitude and direction 
of the spin components, during the merger [bH . Il5l . Il6j |. 
In addition, the radiation of net linear momentum by a 
black-hole binary leads to the recoil of the final remnant 

^eU^UMmMMMMM a in he M 

M, M, SO, 5E 52, 5E 0, HE 5a 51,which can have 
astrophysically important effect s0, [J 5E 5E [IE [IE ■ 

In [2| we introduced the following heuristic model for 
the gravitational recoil of a merging binary. 

V rccoi i{q,di) = v m ei + uj_(cos(£) e\ + sin(£) e 2 ) + W|| e z , 

(1) 
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A = 1.2 x 10 4 kms" 1 HE, B = -0.93 [3J 
(6.9 ± 0.5) x 10 3 kms -1 , on — Si/mf, S t and to, are 
the spin and mass of hole i, q — mi/rri2 is the mass ra- 
tio of the smaller to larger mass hole, 77 = q/(l + q) 2 
is the symmetric mass ratio, the index _L and |j refer to 
perpendicular and parallel to the orbital angular momen- 
tum respectively at the effective moment of the maximum 
generation of the recoil (around merger time), ei,e2 are 
orthogonal unit vectors in the orbital plane, and £ mea- 
sures the angle between the "unequal mass" and "spin" 
contributions to the recoil velocity in the orbital plane. 
The angle O was defined as the angle between the in- 
plane component of A = (mi + 1712) (S^/ms — Si/mi) and 
the infall direction at merger. The form of Eq. (|2aj) was 
proposed in [sE [52|. while the form of Eqs. (l2lt|) and |2c|) 
was proposed in Q based on the post-Newtonian ex- 
pressions in [53|. In Ref Q we determined that K = 
(6.0±0.1) x 10 4 kms -1 , and made the first prediction that 
the maximum possible recoil is ~ 4000 kms -1 for equal- 
mass binaries with anti-parallel spins in the orbital plane 
(in Ref. [4?| . we performed simulations with a measured 
recoil of 3250 kms -1 ). Although £ may in general de- 
pend strongly on the configuration, the results of [37j and 
post-Newtonian [53j calculations show that £ is 90° for 
headon collisions, and the results presented in Ref. [54| 
indicate that £ ~ 145° for a wide range of quasi-circular 
configurations. A simplified version of Eq. |T]) that mod- 
els the magnitude of K- eC oii was independently proposed 
in 139*] , and a simplified form of Eq. (|2bD for the equal- 
mass aligned spin case was proposed in [36j . A more gen- 
eral formula, using only symmetry arguments and fits to 
numerical data, was recently proposed in (55l . [HE] • 

Our heuristic formula |T]) describing the recoil veloc- 
ity of a black-hole binary remnant as a function of the 
parameters of the individual holes has been theoretically 
verified in several ways. In [3J] the cos O dependence was 
established and was confirmed in [43J|for binaries with 
different initial separations. In Ref. [421 ] the decomposi- 
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tion into spin components perpendicular and parallel to 
the orbital plane was verified, and in [46j it was found 
that the quadratic-in-spin corrections to the in-plane re- 
coil velocity are less than 20 kms -1 . Recently, Baker et 
al. [57| measured the recoil for unequal-mass, spinning 
binaries, with spins lying in the initial orbital plane, and 
concluded that the leading order dependence of the out- 
of-plane kick was 0(n 3 ), rather than the 0(r} 2 ) that we 
predicted. In this paper we examine this dependence in 
detail. 

As pointed out in (57J the consequences of an 0{rf) de- 
pendence of the recoil, rather than an 0(rj 2 ) dependence, 
are significant for both the retention of intermediate mass 
black holes (IMBH) in globular clusters and supermas- 
sive black holes in galaxies. It is thus important that we 
understand how the recoil depends on mass ratio. 

In their paper, Baker et al. [57| analyzed configura- 
tions that required fitting two angle parameters (for a 
given value of q) before the maximum recoil could be 
obtained. Their configuration also produced very small 
recoil velocities for smaller values of q. In order to help 
clearly display the dependence, we choose configurations 
that only require fitting one angle parameter in Eq. (|2c[) . 
and have a substantial recoil even for our smallest mass 
ratios. In addition, the Baker et al. [13] runs were con- 
structed such that the spin of the larger black hole is 
proportional to the mass ratio squared times the spin of 
the smaller BH, i.e. a-i = q 2 a±, where q — raxjra^ < 1, 
and cci < 1 is the spin of the smaller black hole. In the 
small mass ratio limit, this leads to a small, very highly 
rotating hole (for any spin-induced recoil to be observ- 
able) orbiting a large, essentially nonrotating, BH. This 
configuration does not match the expected astrophysical 
scenario since large BHs are expected to have high spins. 
Our choice of configuration more closely matches the as- 
trophysical scenario when q — > because in that limit 
the recoil becomes independent of the spin of the smaller 
BH, as is apparent in the factor of a2 — qot\ in Eq. ([2c]) . 

The paper is organized as follows, in Sec. [IT] we review 
the numerical techniques used for the evolution of the 
black-hole binaries and the analysis of the physical quan- 
tities extracted at their horizons, in Sec. IIIII we present 
results and analysis, and in Sec. II VI we present our con- 
clusions. 



II. TECHNIQUES 

To c omp ute initial data, we use the puncture ap- 
proach [58| along with the TwoPunctures [59| thorn. 
In this approach the 3-metric on the initial slice has 
the form jab = (jPbl + «) 4 <5ab, where iJjbl is the Brill- 
Lindquist conformal factor, S a b is the Euclidean met- 
ric, and u is (at least) C 2 on the punctures. The 
Brill-Lindquist conformal factor is given by ipBL = 1 + 
y^_, m^/(2|r— fi\), where n is the total number of 'punc- 
tures', is the mass parameter of puncture i (mf is not 
the horizon mass associated with puncture i), and f$ is 



the coordinate location of puncture i. We evolve these 
black-hole-binary data-sets using the LazEv 16011 imple- 
mentation of the moving puncture approach [lOl [llj ] . In 
our version of the moving puncture approach we replace 
the BSSN [6l|, [62|, [63| conformal exponent </>, which has 
logarithmic singularities at the punctures, with the ini- 
tially C 4 field x = exp(— 40). This new variable, along 
with the other BSSN variables, will remain finite pro- 
vided that one uses a suitable choice for the gauge. An 
alternative approach uses standard finite differencing of 
4> [HI . Recently Marronetti et al. [64| proposed the use of 
W = y/x as an evolution variable. For the runs presented 
here we use centered, eighth-order finite differencing in 
space [65[ and an RK4 time integrator (note that we do 
not upwind the advection terms). 

We use the CACTUS framework @ with the CAR- 
PET [67| mesh refinement driver to provide a 'moving 
boxes' style mesh refinement. In this approach refined 
grids of fixed size are arranged about the coordinate cen- 
ters of both holes. The CARPET code then moves these 
fine grids about the computational domain by following 
the trajectories of the two black holes. 

We obtain accurate, convergent waveforms and horizon 
parameters by evolving this system in conjunction with 
a modified 1+log lapse and a modified Gamma-driver 
shift condition [id . |68|, and an initial lapse ait = 0) = 
2/(1 + ijj% L ). The lapse and shift are evolved with 

{d t ~(3 l d l )a = ~2aK (3a) 
d t f3 a = B a (3b) 
d t B a = 3/4d t f a - oB a . (3c) 

Note that we denote the Gamma-driver parameter by a 
rather than the more typical r\ to avoid confusion with 
the symmetric mass ratio parameter. These gauge con- 
ditions require careful treatment of \, the inverse of the 
three-metric conformal factor, near the puncture in order 
for the system to remain stable (lOl [l2ll20l|. In our tests, 
W showed better behavior at very early times [t < 10M) 
(i.e. did not require any special treatment near the punc- 
tures), but led to evolutions with lower effective resolu- 
tion when compared to x- Interestingly, a mixed evo- 
lution system that evolved W for t < 10M and x f° r 
t > 10M showed inaccuracies similar to the pure W sys- 
tem. At higher resolution W and x agreed with good 
accuracy. In Ref. [6!| it was shown that this choice of 
gauge leads to a strongly hyperbolic evolution system 
provided that the shift does not become too large. 

We use AHFinderdirect [7(3] to locate apparent hori- 
zons. We measure the magnitude of the horizo n sp in us- 
ing the Isolated Horizon algorithm detailed in [71J . This 
algorithm is based on finding an approximate rotational 
Killing vector (i.e. an approximate rotational symmetry) 
on the horizon, and given this approximate Killing vector 
tp a , the spin magnitude is 

S [v] = ^<f AH (<p a R b K ab )d 2 V (4) 



3 



where K a f, is the extrinsic curvature of the 3D-slice, d 2 V 
is the natural volume element intrinsic to the horizon, 
and R a is the outward pointing unit vector normal to 
the horizon on the 3D-slice. We measure the direction of 
the spin by finding the coordinate line joining the poles 
of this Killing vector field using the technique introduced 
in [l6[ • Our algorithm for finding the poles of the Killing 
vector field has an accuracy of ~ 2° (see [16( for details). 
The mass of the horizon is given by the Christodoulou 
formula 

m H = yJm*„ + Sy(*m*„), (5) 

where m m is the irreducible mass. 

We also use an alternative quasi-local measurement of 
the spin and linear momentum of the individual black 
holes in the binary that is based on the coordinate ro- 
tation and translation vectors [45| . In this approach the 
spin components of the horizon are given by 

5 H = ir/ 4>l]R h K ah d 2 v, (6) 

87r J AH 

where = 5£j5 m kr m e l: > k , and r m = x m — x™ is the 
coordinate displacement from the centroid of the hole, 
while the linear momentum is given by 

p ^ = hj AH ^ Rb{Kab ~ K ~i<*)* v > ^ 

where £L = 8\. 

We measure radiated energy, linear momentum, and 
angular momentum, in terms of ^4, using the formulae 
provided in Refs. [72i,l73T|. However, rather than using the 
full tpi we decompose it into £ and m modes and solve 
for the radiated linear momentum, dropping terms with 
I > 5 (only the £ = 2 and £ = 3 modes make significant 
contributions to the recoil). The formulae in Refs. [72|.|73| 
are valid at r = oo. We obtain highly accurate values 
for these quantities by solving for them on spheres of 
finite radius (typically r/M = 50, 60, • • • ,100), fitting 
the results to a polynomial dependence in I = 1/r, and 



extrapolating to I — 0. We perform fits based on a linear 
and quadratic dependence on I, and take the final values 
to be the quadratic extrapolation with the differences 
being the extrapolation error. 



A. Initial Data 

We evolve quasi-circular configurations with a more 
massive spinning black hole, with specific spin a/m H = 
0.8 pointing in the initial orbital plane, and a non- 
spinning smaller BH. The orbital parameters where cho- 
sen using 3PN parameters for quasi-circular orbits with 
orbital period Mlo = 0.05, which provides the puncture 
locations and momenta. We normalize the puncture mass 
parameters so that the total ADM mass is 1M and the 
mass ratio is the specified one. We then modify the con- 
figurations by rotating the initial spin direction (which 
has a small effect on the total ADM mass). We denote 
these configurations by QXXXTHYYY, where XXX de- 
notes the mass ratio (XXX=100 for q — 1, XXX=66 for 
q = 2/3, XXX=50 for q = 1/2, XXX=40 for q = 1/2.5, 
XXX=33 for q = 1/3, XXX=25 for q = 1/4, XXX=17 for 
q = 1/6, and XXX=13 for q = 1/8) and YYY gives the 
angle (in degrees) between the initial spin direction and 
the y-axis. The initial data parameters are summarized 
in Table HI These configuration have several advantages 
when modeling the out-of-plane kick as a function of mass 
ratio. First, one need only fit to two parameters (one an- 
gle and one amplitude) to determine the maximum recoil 
for a given mass ratio; improving the statistical reliability 
of the estimated maximum recoil. Second, the maximum 
recoil velocity is quite large, even for q = 1/6, ensur- 
ing that errors in measuring the maximum recoil veloc- 
ity are not a significant fraction of the recoil itself. In 
addition the functional form used in our non-linear fits, 
Acos(z9 — B), yields a more robust measurement of A 
when compared to fits of Ccos($i — D) + i?cos(z?2 — F) 
(as used in Baker et al.) for small sample sizes (i.e. the 
functional form used here is more amenable to an accu- 
rate fit when the sample size is small). (See Fig. [SJ 



B. Determining the Orbital Plane 

These configuration show significant orbital precession, 
as demonstrated in Fig. [TJ that presents a significant 
challenge when modeling the recoil. Our empirical for- 
mula ([1]) decomposes the recoil in terms of velocities 
parallel and perpendicular to the angular momentum. 
Thus we need an accurate determination of the orbital 
plane near merger (where most of the recoil is gener- 
ated [H, H3|)- We determine an approximate plane by 
choosing points along the plunge trajectory near where 
|r| and |r| have extrema. We then find a rotation, such 



that the late-time approximate orbital plane is rotated 
onto the xy-plane. In order to model the out-of-plane 
recoil, we need to measure it as a function of the orien- 
tation of the spin vector during merger. We do this by 
fixing the remaining freedom in the transformation such 
that the transformed trajectories (for a given sequence 
of configurations) coincide (approximately) near merger. 
We then measure the spin direction at some fixed fidu- 
cial point along the merger trajectory (in practice, at the 
point ro discussed below) and fit the out-of-plane kick to 
the form V = A cos(i9 — B), where ■& is the angle between 
the in-plane component of the spin at this fiducial point 



4 



TABLE I: Initial data parameters for the quasi-circular configurations with a non-spinning smaller mass black hole (labeled 
1), and a larger spinning black hole (labeled 2). The punctures are located at fi = (xi, 0, 0) and f*i = (x2, 0, 0), with momenta 
P — ±(0, P, 0), spins S2 = (S x , S v ,0), mass parameters m p , horizon (Christodoulou) masses m H , and total ADM mass Madm- 
The configuration are denoted by QXXXTHYYY where XXX gives the mass ratio (0.17, 0.25, 0.40, 0.33, 0.50, 0.66, 1. 00) 
and YYY gives the angle in degrees between the initial spin direction and the y-axis. In all cases the initial orbital period is 
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,8072 


Q66THI30 


3 


,93712 


-2 


.61286 


0. 


,128421 


0. 


,38234 





.36743 


-0.22690 


-0.19039 





,40388 





,60582 


0. 


99985 





,8071 


Q66TH2I0 


3 


,93712 


-2 


.61286 


0. 


,128421 


0. 


.38234 





.36743 


0.14810 


-0.25651 





,40389 





,60585 


0. 


99994 





,8070 


Q66TH315 


3 


,93712 


-2 


,61286 


0. 


,128421 


0. 


.38234 





.36743 


0.20944 


0.20944 





,40389 





,60583 


0. 


99988 





,8071 


Q100TH000 


3 


,28027 


-3 


.28027 


0. 


,133568 


0. 


,48338 





.30398 


0.00000 


0.20595 


0, 


,50543 





,50547 


1. 


00001 





,8061 


QI00TH090 


3 


,28027 


-3 


.28027 


0. 


,133568 


0. 


,48338 





.30398 


-0.20595 


0.00000 


0, 


,50541 





,50537 


0. 


99965 





,8065 


QI00TH130 


3 


,28027 


-3 


,28027 


0. 


.133568 


0. 


.48338 





.30398 


-0.15777 


-0.13238 





,50542 





,50541 


0. 


99980 





,8063 


QI00TH210 


3 


,28027 


-3 


,28027 


0. 


.133568 


0. 


,48338 





.30398 


0.10297 


-0.17836 





,50543 





,50544 


0. 


99992 





,8062 


QI00TH3I5 


3 


,28027 


-3 


.28027 


0. 


,133568 


0. 


,48338 





.30398 


0.14563 


0.14563 





,50542 





,50542 


0. 


99983 





,8063 



for a given configuration with the in-plane spin for the 
corresponding QXXXTH000 configuration. 

We determine the orbital plane during merger in the 
following way. We plot r = \r\ — r*2| as a function of 
time and determine the location of the maximum of |f|. 



We then choose two points on either side of the maxi- 
mum with similar values of |f| and a third point close 
to the maximum. We denote these three points of the 
trajectory with r+, r*o, an d r_. Wc then choose a rota- 
tion such that the vector f a = r + — f lies on the new 
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THOOO' 

TH210' 




FIG. 1: The trajectory fi - r 2 for the Q25TH000 and 
Q25TH210 configurations. Note the significant precession and 
the lack of alignment between the late-time orbital planes. 

y-axis and that the normal to the plane determined by 
r a and ft, = r _ — r*o lies along the new z-axis. With these 
choices we uniquely determine a rotational transforma- 
tion from the coordinates used by the code to coordinates 
where the plunging orbital plane coincides with the xy 
plane. This transformation also rotates the trajectories 
in the appropriate way such that trajectories with similar 
late-time dynamics will overlap during merger (this is a 
generalization of the procedure given in Ref. [471. l54j'). 

In Fig. [1] we show the orbital trajectory difference 
fi - f 2 for the Q25TH000 and Q25TH210 configurations. 
Note the significant precession and that the orbital planes 
near merger do not coincide. Fig. [2] shows that, after ro- 
tating the plane, as described above, the merger trajec- 
tories coincide. 

In Figs. [3] and 2] we show the xy projections of the tra- 
jectories for the Q25THYYY configurations before and 
after rotating the merger orbital plane. Note that prior to 
this transformation, there is no rotation in the xy plane 
that will make the late-time trajectories overlap and the 
reasonable overlap of the late-time trajectories for the 
transformed case. 



III. RESULTS 

We used the following grid configurations for the runs. 
For the q = 1/2, q — 1/2.5, and q — 1/3 runs we placed 
the outer boundaries at 1664M with a coarsest resolution 
of h = 25. 6M, we used 12 levels of refinement around the 
smaller BH, with finest resolution of h = M/80, and 11 
levels of refinement around the larger BH. For the q = 1/4 
runs we added an additional level of refinement around 
the smaller BH, for q = 1/6 and q = 1/8 we added two 
additional levels (for a total of 14), and for the q = 2/3 
runs we added an additional level of refinement about 
the larger BH. For the q = 1 runs we used 11 levels of 



THOOO 

TH210 




FIG. 2: The trajectory fi - r 2 for the Q25TH000 and 
Q25TH210 configurations after rotating the system. Note the 
good agreement in the late time trajectories and that these 
trajectories now lie on the xy plane. The solid points are the 
locations of r+, r*o, and r_. Note that these points agree for 
both curves. 




■ TH315 

-8 I ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' 

-8 -6 -4 -2 2 4 6 8 

(x, - x 2 )/M 

FIG. 3: The xy projection for the untransformed trajectory 
fi — ?2 for the Q25THYYY configurations. The trajectories 
have been rotated (in the xy plane). Note that there are no 
rotations in the xy plane that will make the late-time trajec- 
tories overlap. The filled circles are the locations of r+, fo, 
and r_. 



refinement around the non-spinning hole and 12 levels 
of refinement around the spinning hole. In the q = 1/4 
case, we found that using eighth-order accuracy was crit- 
ical, as incorrect dynamics (outspiral rather than inspiral 
when evolving with x, and plunges rather than inspirals 
when evolving with W) resulted when using fourth or- 
der methods. The dynamics obtained from the W and 
X systems agreed when using eighth-order methods. In 
all cases we used a = 3 in the Gamma-driver shift con- 
dition. In addition we ran all q — 1/4 configuration with 
a higher resolution (by a factor of 1.2) for all levels (See 



FIG. 4: The xy projection for the transformed trajectory 
fi — ?2 for the Q25THYYY configurations. No additional 
rotations have been applied. Note the good agreement in the 
late time trajectories. The filled circles are the locations of 
r+ , fo , and r_ . 



In Table |TT] we give the radiated energy and angular 
momenta for all configurations, as well as the untrans- 
formcd recoil velocities. 



In order to measure the dependence of our recoil cal- 
culations on the choice of r + , ro, r_ (the three points 
in the trajectory that define the late-time orbital plane) 
we make two choices for these quantities a fiducial 
choice (r + ,ro,r_)/M = (2.0, 1.0, 0.5) and a second choice 
(r+, ro, r~)/M = 2.2,1.2,0.7) based on the point in the 
orbital trajectory when \r(r)\ is a maximum. For our 
runs we found that the maximum in \r(t)\ occurs at 
r = (1.2 ±0.1). We also use an alternative approach, 
discussed in Sec. lIII Cl where we choose (r + , r , r_) to be 
the extrema of f(t) and f(t) during the plunge phase. In 
this latter approach, the values of (r+,ro,r_) vary with 
each configuration and generally produce poorer results 
than for the fixed choice of (r+ , ro , r_) . In Table IIIII we 
show how these two choices affect the calculation of 



and v± , the spin component in and out of the plane, and 
the angle between the in-plane spins with the correspond- 
ing spin for the TH000 configurations (i9). Note that we 
report the magnitude of the in-plane (i.e. v±) recoil and 
spin aj_, because modifying (r + ,ro,r_) introduces addi- 
tional rotations within the orbital plane; making a direct 
comparison of x and y components meaningless. Note 
that there is scatter in both the in-plane and out-plane 
components of the spin. The scatter in the out-of-plane 
component appears to be relatively significant due to its 
small average value. The source of this scatter may sim- 
ply be due to errors in estimating the true direction of 
the orbital plane (since the spin in the plane is large). We 
take the average value of the in-plane spin when fitting 
to Eq. flU below. 



Our empirical formula |T]) predicts that vn will scale as 
cos($ — t?o)- For each set of configurations with a given 
mass ratio, we perform a non-linear least squares fit of 
Vii to the form 

W|| = Acos(tf--B) (8) 



and solve for A and B. The results of these fits are sum- 
marized in Table ITVl In Fig.[5]we plot the individual data 
points for each mass ratio and the best fit function. Note 
that the effect of refining the grid is to slightly increase 
the magnitude of the out-of-plane recoil. 



With the above results we can now fit the maximum 
out-of-plane recoil (uii) versus q. Our empirical for- 
mula |T]) predicts that this maximum will have the form 



v\\ = Kri 2 /(\ + q)(ai - qa 2 ), where r\ = q/(l + q) 2 . 
For all of our configurations a 2 ~ 0, and hence Vu = 
aiKrj 2 /(l + q). Interestingly, this form has a maximum 
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TABLE II: The radiated energy and angular momentum, as well as the recoil velocities for each configuration. Note that some 
of the error estimates, which are based on the differences between a linear and quadratic extrapolation in I — 1/r, are very 
small. This indicates that the differences between the extrapolation can underestimate the true error. All quantities are given 
in the coordinate system used by the code (i.e. the untransformed system). The Q25THXXX quantities given on the bottom 
of the table are for the higher resolution runs. See Sec. IIII Al for an explanation of the differences in the radiated quantities. 



Config 


WOE/M 


100 J x /M 2 


100 Jy/M 2 


100 J z /M 2 


V x 


Vy 


v z 


Q13TH000 


0.623 


± 


0.018 


2.383 ±0.031 


1.987 ±0.050 


2.002 


± 


0.035 


66.5 ±2.3 


26.9 ±5.0 


101.0 ±6.8 


Q13TH090 


0.6444 


:=t 


: 0.005 


-2.182 ±0.008 


2.533 ±0.017 


2.142 


± 


0.184 


-14.2 ± 6.5 


70.7 ±3.3 


92.7 ±6.3 


Q13TH130 


0.671 


± 


0.013 


-3.438 ± 0.034 


0.509 ± 0.008 


2.112 


± 


0.012 


-393.0 ±11.9 


-167.4 ± 15.3 


43.2 ±0.1 


Q13TH210 


0.637 


± 


0.007 


-1.072 ±0.031 


-2.980 ±0.038 


2.199 


± 


0.140 


38.1 ±3.9 


86.3 ±2.8 


-97.1 ±7.5 


Q13TH315 


0.689 


± 


0.011 


3.484 ± 0.041 


-0.161 ±0.009 


2.162 


± 


0.049 


-382.0 ± 12.3 


-216.2 ±17.3 


-27.2 ±0.2 



Q17TH000 


l.i 


032: 


±i 




3343 


3.644 ± 0.027 


3.990 


±0.094 


3.70 ± 0.094 


-99.5 ±10.8 


-624.5 ±27.6 


-137.9 ±6.8 


Q17TH090 





.980 


± 





.005 


-4.133 ±0.014 


3.903 


± 0.078 


3.66 ± 0.246 


323.1 ±5.2 


-105.5 ±1.8 


-146.1 ± 10.0 


Q17TH130 


1 


.049 


± 





.027 


-5.592 ±0.059 


0.004 


±0.014 


3.484 ± 0.045 


371.8 ± 10.4 


343.2 ± 16.4 


-42.3 ± 0.3 


Q17TH210 


1 


.027 


± 





.032 


-1.162 ±0.011 


-5.351 ±0.100 


3.70 ± 0.077 


228.6 ± 21.5 


-579.2 ± 14.8 


140.8 ± 5.4 


Q17TH315 


1 


.032 


± 





.012 


5.401 ± 0.044 


0.352 


±0.035 


3.687 ±0.222 


175.0 ±8.3 


270.0 ± 10.6 


-30.3 ± 3.5 



Q25TH000 


1 


.621 


± 





.020 


4.543 ±0.135 


5.521 ±0.044 


7.571 


± 


0.115 


145.1 ± 15.6 


855.3 ±21.9 


535.3 ±9.9 


Q25TH090 


1 


.642 


± 





.016 


-5.895 ±0.076 


4.412 ±0.039 


7.868 


± 


0.233 


-869.9 ± 19.8 


80.5 ± 14.2 


446.3 ± 7.4 


Q25TH130 


1 


.530 


± 


0. 


.012 


-6.879 ± 0.040 


-0.096 ± 0.005 


7.400 


± 


0.209 


-64.7 ±9.1 


-141.3 ±3.5 


-131.7 ±2.2 


Q25TH210 


1 


.638 


± 


0. 


.018 


-1.045 ±0.106 


-7.132 ±0.021 


7.776 


± 


0.062 


-335.9 ± 27.8 


819.1 ± 17.2 


-487.8 ± 8.5 


Q25TH315 


1 


.512 


± 





.014 


6.774 ± 0.063 


0.583 ±0.030 


7.375 


± 


0.188 


73.4 ±7.1 


-45.9 ± 0.4 


223.4 ±3.9 



Q33TH000 


2 


.216 


± 


0. 


.023 


4.731 ± 0.073 


7.171 ±0.087 


11.516 


± 


0.098 


47.5 ± 19.1 


983.3 ± 19.6 


827.3 ±8.8 


Q33TH090 


2 


.184 


± 


0. 


.018 


-7.569 ± 0.026 


4.653 ± 0.040 


11.548 


± 


0.082 


-829.1 ± 21.2 


-40.4 ± 15.5 


576.6 ± 8.0 


Q33TH130 


2 


.050 


± 





.016 


-8.485 ±0.017 


-0.840 ± 0.053 


11.108 


± 


0.089 


171.0 ±4.8 


8.7 ±0.2 


-374.3 ± 0.3 


Q33TH210 


2 


.200 


± 





.020 


-0.433 ± 0.081 


-8.597 ±0.062 


11.515 


± 


0.032 


-453.0 ± 26.2 


793.0 ± 12.0 


-697.3 ± 6.7 


Q33TH315 


2 


.054 


± 





.017 


8.335 ± 0.022 


1.495 ±0.010 


11.072 


± 


0.002 


306.9 ± 5.6 


155.7 ±5.4 


520.5 ±0.7 



Q40TH000 


2 


.613 


± 





.021 


4.604 ± 0.114 


8.337 ±0.024 


14.791 


± 


0.074 


19.6 ± 19.4 


985.0 ± 15.0 


1132.5 ±4.8 


Q40TH090 


2 


.607 


± 





.017 


-8.761 ±0.001 


4.729 ±0.133 


14.790 


± 


0.147 


-929.8 ± 18.2 


-83.4 ± 16.3 


939.3 ± 6.4 


Q40TH130 


2 


.452 


± 


0. 


.012 


-9.332 ± 0.033 


-1.734 ±0.082 


14.123 


± 


0.130 


-12.3 ± 1.7 


-168.0 ±5.5 


-137.0 ±3.6 


Q40TH210 


2 


.621 


± 


0. 


.018 


0.157 ±0.047 


-9.604 ±0.013 


14.815 


± 


0.072 


-520.1 ± 24.9 


845.9 ± 10.4 


-1068.0 ±5.8 


Q40TH315 


2 


.441 


± 


0. 


.014 


9.165 ±0.068 


2.464 ±0.083 


13.933 


± 


0.039 


120.6 ± 1.3 


-39.4 ±0.6 


328.5 ±2.9 



Q50TH000 


2. 


.814 


± 





.013 


4.089 ± 0.0848 


7.762 ±0.092 


16.672 ± 0.023 


-109.8 ±3.1 


159.5 ± 7.5 


-77.3 ± 8.8 


Q50TH090 


2. 


.865 


± 


0. 


.018 


-8.094 ± 0.0421 


4.566 ±0.198 


17.499 ± 0.336 


623.4 ± 1.9 


-75.4 ± 13.3 


-1118.0 ±7.8 


Q50TH130 


2. 


.991 


± 





.017 


-9.082 ± 0.1534 


-2.150 ±0.025 


17.597 ±0.114 


688.2 ± 6.8 


611.5 ± 20.9 


-1330.0 ± 1.0 


Q50TH210 


2. 


.805 


± 





.014 


0.087 ± 0.0538 


-8.922 ± 0.098 


16.83 ±0.089 


9.5 ±4.9 


-221.3 ±5.7 


590.5 ± 9.9 


Q50TH315 


2. 


985 


± 





.018 


8.788 ±0.1391 


2.945 ± 0.0222 


17.43 ±0.0177 


600.4 ± 8.4 


665.1 ±20.4 


1253.4 ±2.4 



Q66TH000 


3 


.277 


± 





.012 


3.404 ±0.1164 


7.633 ± 0.0533 


20. 


.177 


± 





.193 


-69.7 


± 


0.9 


-294.6 ±2.2 


-848.5 ±9.2 


Q66TH090 


3 


.419 


± 





.013 


-8.028 ± 0.0443 


3.166 ±0.2086 


21. 


293 


± 





.133 


749.2 


± 


6.2 


17.3 ± 11.4 


-1565.8 ±11.9 


Q66TH130 


3 


.344 


± 





.012 


-7.896 ± 0.1482 


-2.756 ±0.2232 


20. 


.736 


± 





.060 


356.3 


± 


1.0 


425.9 ± 11.8 


-940.4 ± 2.3 


Q66TH210 


3 


.336 


± 


0. 


.013 


0.927 ±0.1165 


-8.427 ±0.0148 


20. 


701 


± 





.048 


221.8 


± 


7.5 


-487.5 ± 2.9 


1265.4 ± 11.0 


Q66TH315 


3 


.321 


± 


0. 


.011 


7.628 ±0.1358 


3.368 ±0.1800 


20. 


.651 


± 





.087 


259.2 


± 


3.6 


404.5 ± 11.2 


750.7 ± 3.2 



Q100TH000 


3.580 


± 


0.006 


2.036 ±0.1025 


6.224 ± 0.0869 


23.232 ± 0.038 


-10.2 ± 2.6 


-447.0 ± 2.6 


-1354.5 ± 5.2 


Q100TH090 


3.538 


± 


0.004 


-6.452 ± 0.2779 


2.110 ±0.2082 


23.541 ± 0.332 


344.2 ± 11.2 


31.1 ±3.3 


-962.9 ±4.6 


Q100TH130 


3.461 


± 


0.001 


-6.393 ± 0.3509 


-2.042 ±0.3155 


22.749 ± 0.003 


-46.7 ± 14.3 


6.8 ±3.8 


245.1 ± 12.6 


Q100TH210 


3.606 


± 


0.011 


1.268 ± 0.0057 


-6.940 ± 0.5645 


23.382 ± 0.072 


241.4 ±0.6 


-418.6 ±0.3 


1406.9 ±2.8 


Q100TH315 


3.466 


± 


0.001 


6.268 ± 0.4289 


2.608 ±0.2515 


22.748 ± 0.035 


-77.6 ± 12.7 


-34.8 ±3.4 


-383.3 ± 12.9 


Q25TH000HR 


1.605 


± 


0.022 


4.623 ± 0.025 


5.699 ±0.014 


7.503 ±0.124 


110.2 ±3.1 


539.2 ±0.2 


471.9 ±4.6 


Q25TH090HR 


1.632 


± 


0.022 


-6.164 ±0.050 


4.759 ±0.049 


7.668 ± 0.008 


-698.5 ± 1.7 


104.1 ±6.2 


504.7 ±4.5 


Q25TH130HR 


1.736 


± 


0.021 


-7.811 ±0.057 


-0.682 ±0.012 


7.793 ± 0.012 


-683.3 ± 2.4 


-549.1 ± 13.4 


400.6 ± 1.6 


Q25TH210HR 


1.629 


± 


0.022 


-1.129 ±0.002 


-7.430 ± 0.038 


7.622 ±0.134 


-245.6 ±8.2 


656.9 ±0.9 


-512.7 ±4.4 


Q25TH315HR 


1.715 


± 


0.020 


7.618 ± 0.043 


1.323 ±0.023 


7.692 ± 0.005 


-539.0 ± 1.7 


-549.2 ± 11.6 


-300.4 ±0.7 
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TABLE III: The transformed spin (at r = ro), recoil velocities in kms -1 , and angles between the in-plane QXXXTH000 spins 
and the other QXXXTHYYY configurations in degrees. A '1' denotes a quantity measured using (r+, ro, r_)/M = (2.0, 1.0, 0.5), 
while a '2' denotes a quantity measured using (r+, ro, r_)/A/ = (2.2,1.2,0.7). The differences between the values given for 
(2.0, 1.0,0.5) and (2.2, 1.2,0.7) are indicative of the errors. Note that in the transformed system L z and S z are both negative 
(i.e. there is some partial spin/orbit alignment. We denote quantities in the (transformed) orbital plane with a _L subscript. 
The Q25THXXX quantities given on the bottom of the table are for the higher resolution runs. See Sec. MI Al for an explanation 
of the differences. The time derivative of the orbital separation rj is a maximum at r ~ 1.2. 



Config 


v±(i) 


*Tl(l) 


a±/m B (l) 


a|,/m H (l) 


tf(l) 


V ± {2) 


Ml (2) 


a ± /m H (2) 


a l|. 


lm H (2) 


m 


Q13TH000 


57.1362 


-109.88 


0.768518 


-0.237098 


0. 


49.953 


-113.327 


0.751038 


-0. 


.291919 


0. 


Q13TH090 


64.941 


-97.8627 


0.773066 


-0.22036 


3.63707 


58.4083 


-101.897 


0.756446 


-0 


.277112 


3.96937 


Q13TH130 


220.34 


-368.446 


0.727118 


-0.332954 


-78.283 


249.903 


-349.071 


0.761071 


-0. 


.249773 


-77.2282 


Ol 3TH91 


42.1564 


1 98 67Q 


763909 


-0.24517 


177.413 


31 Q5Q4 


131 589 


746579 


-0. 


9Q5Q65 


177.276 


Q13TH315 


239.084 


369.106 


0.742796 


-0.295281 


87.2139 


280.513 


338.693 


0.772431 


-0 


.212227 


89.2203 


Q17TH000 


309.622 


568.413 


0.751456 


-0.270601 


0. 


370.742 


530.574 


0.776806 


-0. 


.193091 


0. 


Q17TH090 


134.116 


344.786 


0.735493 


-0.317919 


67.034 


123.344 


348.784 


0.738562 


-0. 


.314279 


63.9956 


Q17TH130 


351.795 


366.112 


0.78253 


-0.15129 


-36.082 


393.985 


320.27 


0.786039 


-0. 


.145002 


-35.7209 


Ol 7TH91 D 






D 74fS0Q' : i 


n 97QD1 1 


1 7fi ^1 1 


Qc;7 ORG 




n 77^07^ 


n 
u 


1 QQ9QQ 


1 7fi d<\Q 


Q17TH315 


283. 


-156.073 


0.784854 


-0.145801 


123.407 


295.546 


-130.769 


0.778357 


-0 


.186578 


123.21 


Q25TH000 


373.905 


-948.305 


0.757314 


-0.240225 


0. 


459.51 


-909.911 


0.77861 


-0 


.167965 


0. 


Q25TH090 


452.721 


-870.357 


0.779284 


-0.154404 


-21.2857 


544.717 


-815.942 


0.788089 


-0. 


.114385 


-21.0555 


Q25TH130 


197.85 


48.3342 


0.770647 


-0.202858 


-84.2701 


195.719 


56.3483 


0.759915 


-0 


.247894 


-85.4253 


Q25TH210 


426.136 


916.571 


0.77451 


-0.184646 


165.452 


521.338 


865.968 


0.787368 


-0 


.131325 


165.72 


Q25TH315 


124.649 


-204.556 


0.759673 


-0.241429 


85.9471 


115.438 


-209.892 


0.750901 


-0. 


.274232 


84.2327 


Q33TH000 


477.923 


-1193.82 


0.787385 


-0.136406 


0. 


585.012 


-1145.16 


0.794072 


-0 


.101676 


0. 


Q33TH090 


506.036 


-874.874 


0.794658 


-0.0857209 


-23.2403 


573.118 


-832.474 


0.794033 


-0 


.100666 


-23.2941 


Q33TH130 


121.993 


393.11 


0.754769 


-0.264995 


-83.3346 


106.146 


397.682 


0.753714 


-0. 


.275061 


-85.3148 


Q33TH210 


510.834 


1029.22 


0.793461 


-0.100511 


165.757 


597.049 


981.726 


0.795102 


-O.i 


3964032 


165.86 


Q33TH315 


142.073 


-607.563 


0.74435 


-0.292739 


85.1084 


124.37 


-611.433 


0.751109 


-0 


.280681 


82.8364 


Q40TH000 


395.545 


-1447.96 


0.782576 


-0.177756 


0. 


501.737 


-1414.67 


0.792558 


-0 


.121423 


0. 


Q40TH090 


506.998 


-1223.4 


0.797605 


-0.0811617 


-29.708 


599.942 


-1180.61 


0.798469 


-0. 


.074772 


-29.6255 


Q40TH130 


199.308 


85.9742 


0.773172 


-0.20708 


-88.1121 


196.752 


91.6738 


0.767006 


-0 


.237119 


-89.0041 


Q40TH210 


468.966 


1380.9 


0.794067 


-0.117055 


163.468 


578.675 


1338.63 


0.797642 


-O.i 


3855529 


163.572 


Q40TH315 


138.257 


-323.823 


0.762353 


-0.243806 


82.1551 


127.284 


-328.291 


0.760328 


-0 


.257681 


80.7494 



Q50TH000 


206.386 


29.5953 


0.775871 


-0.197482 


0. 


205.713 


33.9596 


0.771006 


-0.223791 


0. 


Q50TH090 


211.758 


1264.66 


0.760737 


-0.256772 


-53.0652 


228.683 


1261.71 


0.77738 


-0.20449 


-53.2402 


Q50TH130 


426.349 


1560.33 


0.793283 


-0.105245 


-109.529 


530.071 


1528.21 


0.800072 


-0.0682939 


-108.518 


Q50TH210 


117.707 


-619.58 


0.759212 


-0.259 


158.194 


100.605 


-622.586 


0.763761 


-0.251091 


157.494 


Q50TH315 


454.34 


-1472.23 


0.799245 


-0.0784754 


62.2872 


549.566 


-1439.4 


0.799962 


-0.065415 


63.1927 


Q66TH000 


118.274 


893.123 


0.763607 


-0.251388 


0. 


101.713 


895.16 


0.772394 


-0.226846 


0. 


Q66TH090 


268.896 


1714.94 


0.7931 


-0.116783 


-62.1364 


354.418 


1699.33 


0.802482 


-0.066631 


-60.9979 


Q66TH130 


336.909 


1038.86 


0.80055 


-0.0605261 


-111.577 


377.096 


1024.96 


0.797946 


-0.0868091 


-110.371 


Q66TH210 


163.971 


-1364.27 


0.770063 


-0.232589 


157.52 


176.949 


-1362.65 


0.783732 


-0.184392 


157.808 


Q66TH315 


311.388 


-835.099 


0.798347 


-0.0778683 


60.8859 


339.961 


-823.881 


0.795121 


-0.109404 


62.06 


Q100TH000 


88.5695 


1423.66 


0.79471 


-0.158448 


0. 


108.207 


1422.3 


0.800356 


-0.116206 


0. 


Q100TH090 


134.476 


1014.12 


0.806191 


-0.0286805 


-65.6896 


167.057 


1009.26 


0.803238 


-0.0439484 


-65.4866 


Q100TH130 


52.4647 


-244.074 


0.785337 


-0.163234 


-116.365 


42.6164 


-245.985 


0.785246 


-0.170854 


-116.398 


Q100TH210 


95.9383 


-1484.52 


0.805662 


-0.0736974 


151.702 


151.312 


-1479.91 


0.804525 


-0.0475139 


151.917 


Q100TH315 


55.0809 


388.766 


0.78288 


-0.178353 


58.0786 


40.6688 


390.537 


0.78611 


-0.178897 


57.8702 


Q25TH000HR 


238.923 


-684.464 


0.727703 


-0.328834 





235.693 


-685.583 


0.744828 


-0.292768 





Q25TH090HR 


302.537 


-813.638 


0.731821 


-0.318884 


-12.0909 


317.701 


-807.838 


0.754461 


-0.263398 


-11.4989 


Q25TH130HR 


514.081 


-815.226 


0.787921 


-0.127218 


-83.1554 


596.917 


-756.679 


0.791774 


-0.114032 


-80.3582 


Q25TH210HR 


299.288 


815.542 


0.730048 


-0.325618 


167.319 


315.17 


809.537 


0.753828 


-0.269311 


167.829 


Q25TH315HR 


501.346 


656.57 


0.79009 


-0.113807 


85.1511 


560.613 


606.749 


0.789798 


-0.127652 


87.8788 




500 
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TABLE IV: Fit parameters for wy = Acos((i9 - B)tt/180) 
for each q. A is in units of kms -1 and B is in degrees. A 
(1) denotes values obtained from the (r+,ro,r—) = (2, 1,0.5) 
transformation, while a (2) denotes values obtained from the 
(r+,ro,r_) = (2.2,1.2,0.7) transformation. The last row 
shows A and B when the grid is refined by a factor of 1.2 

g A(l) B(l) A(2) B(2) 

1/8 381.2 ± 7.1 107.33 ±0.93 355.7 ±4.4 109.35 ± 0.63 
1/6 577 ±11 15.51 ± 1.18 538 ± 12 17.2 ±1.4 

1/4 981 ± 11 188.03 ±0.70 926.3 ± 5.1 187.61 ± 0.33 

1/3 1297.8 ± 9.2 203.75 ± 0.37 1231.7 ± 9.5 203.22 ± 0.41 

1/2.5 1472.4 ± 9.6 184.88 ± 0.39 1424.1 ± 5.1 184.35 ± 0.21 

1 /2 1640 ± 23 269.96 ± 0.77 1603 ± 15 270.32 ± 0.54 

2/3 1709.1 ± 7.7 301.26 ± 0.26 1686.0 ± 8.6 246.04 ± 0.30 

1 1508.1 ± 8.7 342.69 ± 0.32 1502.5 ± 6.7 342.76 ± 0.25 

TJl 994.6 ±6.4 133.03 ±0.43 942.5 ± 4.6 137.09 ± 0.35 



at q = 2/3. Baker et al. (5?J propose that Eq. ([I]) should 
be modified to v\\ = 4Krf/(l + q)(ai — qa2) (which has a 
maximum at q — 3/4 for fixed ot\ and a^ — 0). In order 
to discriminate between these two possibilities we fit our 
data to the form 

v\\/ a ^G{^) H /{!&{! + q)), (9) 

where we set G — K and solve for H, as well as al- 
lowing both G and H to vary. Here a = — qa^\- 
It is important to note that our fits use eight values 
of q to obtain one or two parameters. In Figs. [6l [3 
and [5] we show results from these fits for the choices 
(r+,r ,r_) = (2,1,0.5) and (r+,r ,r_) = (2.2,1.2,0.7). 
In the figures we plot the predicted recoil velocity us- 
ing both our formula and the one proposed by Baker 
et al. assuming an uncertainty in K twice as large as 
that given in From the plots we can see that the 
recoil velocities agree with our empirical formula much 
better than with the Baker et al. modification. The best 
fit functions have the form Vn/a = Kr/ 1 ' 9 062 /(l + q) 
(with a root mean square error in the predicted v/a of 
68 kms -1 ) and v\\/a = Kt] 2 036±00A6 /(1 + q) (with a root 
mean square error in v/a of 46 kms -1 ) for the choices 
(r + ,r ,r_) = (2,1,0.5) and (r+,r ,r_) = (2.2,1.2,0.7) 
respectively, where K was set to K = 6.0 x 10 4 , and 
V\\/a = (62965 ± 746) ?7 2 - 027±0 - 048 /(l + q) (with a root 
mean square error in the predicted v/a of 35 kms -1 ) 
andw||/a = (62024±555)ry 2127±0 037 /(l±g) (with a root 
mean square error in the predicted v/a of 26 kms -1 ) re- 
spectively when both G and H are varied. The value of 
the constant K = Ggt determined in these latter two fits 
is in reasonable agreement with our previous measure- 
ment of K = (6.0 ± 0.1) x 10 4 . 

A. Finite difference errors 

The effect of finite difference errors are most significant 
for the smaller q runs. Here the main effect regarding our 



FIG. 6: Fit of out-of-plane recoil (Vj|/a) versus mass ratio q 
for the (r+,ro,r_) = (2,1,0.5) choice of orbital plane. The 
top shaded region shows our predicted value for K = (6.0 ± 
0.2) x 10 4 , the lower region shows the prediction based on 
the modification proposed by Baker et al., the lower (green) 
curve shows a fit to rf 1 , while the upper (blue) curve has 
a simultaneous fit to rf 1 and K. Note that the measured 
values of Vy/a, while slightly overshooting the predictions, 
agree much better with our 0(rj 2 ) form. The dots are the 
locations of the data points. 

V/tr 




FIG. 7: Fit of out-of-plane recoil (Vy/a) versus mass ratio 
q for the (r+,ro,r_) = (2.2,1.2,0.7) choice of orbital plane. 
The top shaded region shows our predicted value for K — 
(6.0 ±0.2) x 10 4 , the lower region shows the prediction based 
on the modification proposed by Baker et al., the lower (green) 
curve shows a fit to rf 1 , while the upper (blue) curve has 
a simultaneous fit to rf 1 and K. Note that the measured 
values of V\\/a, while slightly overshooting the predictions, 
agree much better with our 0(r] 2 ) form. The dots are the 
locations of the data points. 

V/a 




calculation is that the $ is rotated with respect to its an- 
alytic value (See e.g. [EH]). This means that, although 
there may be significant errors in both the magnitude 
and direction of the recoil in a particular run, the max- 
imum recoil (as determined by the fit (|8|)) for a given q 
is relatively insensitive to this error. This can be seen 
in the dependence of A in Eq. ([H]) on resolution for the 
q = 1/4 case in Table IIVI This relatively small differ- 
ence is smaller than the mean error in the fits to the 
form ([n]). If we perform the fit to form ^ using the 
high- resolution q = 1/4 results, then the fitting param- 
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FIG. 8: The out-of-plane recoil (Vji/a) versus mass ratio q for 
both choices of (r+,ro,r_) as well as the predictions assum- 
ing a leading if and r/ 3 dependence. The solid lines show the 
prediction assuming a leading -q 2 dependence and an error in 
K of 0.2 x 10 4 kms -1 , while the dotted lines show the pre- 
diction assuming a leading rf dependence and an error in K 
of 0.2 x 10 4 kms" 1 . The ((r+,r ,r_) = (2,1,0.5) data points 
lie above the (2.2, 1.2.0.7) data points. The (2.2, 1.2.0.7) data 
points lie closer to the predicted values, which is consistent 
with the observation that ro = 1.2 is a better approximation 
to the location of the maximum of \r\ for all configurations 
than |r = 1|. The differences in these data points (for a given 
q) is indicative of the true error in calculating the out-of-plane 
recoil. 

2150 



1650 



1150 



650 



150 




Kri7(1 + q) 
4Kti /(1+q) 



0.11 



0.31 



0.51 



0.71 



0.91 



eters change from (G = 62965 ± 746, H = 2.03 ± 0.04) 
to (G = 62908 ± 805, H = 2.01 ± 0.05) and from (G = 
62023 ± 555, H = 2.13 ± .04) to (G = 61974 ± 584, H = 
2.11 ±0.04) for the choices (r+,r ,r_) = (2,1,0.5) and 
(r + ,ro,r_) = (2.2,1.2,0.7) respectively (i.e. the change 
is not statistically significant). 

Note that although we expect the q = 1/6 and q = 
1/8 results to be less accurate than those of the more 
modest mass ratio cases, including these simulations does 
not change the observed 77 dependence. This strongly 
suggests that the 77 dependence is robust against finite 
difference errors. 

We note that very recent work in PN prediction of 
the recoil [74[ and perturbative calculations (a work in 
progress by the authors) show that there is an 77 s cor- 
rection to the leading ?7 2 dependence of the out-of-plane 
recoil. An attempt to fit our data to the functional form 
V/a — (Arj 2 + Bff)/{\ + q) does not yield an accurate 
evaluation of B because the functional form of the recoil 
is insensitive to B at the level of accuracy we obtained. 
That is to say, even if B is large, for example B ~ 39000 
and A ~ 52000 is correspondingly smaller (such that 
Bif + Atj 2 ~ 60000/16 when 77 = 1/4) the maximum 
change in the predicted recoil is of order 10 kms -1 over 
the entire range of 77. We note that both the PN and 
perturbative predictions include a leading order rj 2 de- 
pendence for all configurations including those of Baker 
et al. [H3|. 



B. The In-Plane Recoil 

The empirical formula |T]) provides an accurate pre- 
diction for the large out-of-plane recoil as a function of 
mass ratio (and most importantly, an accurate predic- 
tion of the magnitude of the recoil). Interestingly, the 
in-plane recoil seems to be larger than initially expected. 
Examining Table IIII1 we see that the magnitude of the 
in-plane recoil v± is larger than the predictions of Eq. ([1]) 
(the maximum predicted in-plane recoil velocity for the 
unequal mass cases is about 100 - 200 kms -1 and is dom- 
inated by the unequal-mass, rather than spin, component 
of the recoil; while for the equal-mass case, the in-plane 
recoil is ~ 50 kms -1 ). The excess in-plane recoil can be 
understood in terms of a newly discovered higher-order 
effect (here demonstrated for the first time) that is due 
to the precession of the orbital plane. To understand this 
effect we need to consider how the net recoil is generated. 
As the binary merges, asymmetrical radiation leads to an 
oscillation in the momentum of the center of mass. In the 
more symmetrical 'superkick' configurations, with equal- 
mass and equal- and-opposite spins (with spins in the or- 
bital plane), this leads to the center-of-mass alternately 
recoiling up and then down in an oscillatory manner. The 
total recoil is then determined by where in this cycle the 
plunge occurs. Here this upward/downward oscillation 
is superimposed upon a precession of the orbital plane 
itself. Thus, for example, a strong 'upward' recoil (i.e. a 
recoil along the orbital angular momentum direction) at 
one point during the merger will not be exactly in the op- 
posite direction as the following 'downward' recoil. This 
has two effects, the two recoils can never cancel exactly, 
and a recoil that initially was entirely out-of-plane actu- 
ally has an in-plane component with respect to the final 
plunge orbital plane. Thus both the 'upward' and follow- 
ing 'downward' recoil both contribute, in general, to the 
in-plane recoil. To put it another way, in the precessing 
case, there is a correction to the in-plane-recoil which is a 
function of the history of the instantaneous out-of-plane 
recoil and the extent that the orbital plane precesses 
when the instantaneous recoil is large. However, the re- 
coil is only large near merger, and the orbital plane does 
not precess to a high degree during this short time pe- 
riod (making the correction a higher-order effect). Thus 
the effect, while significant for the much smaller in-plane 
recoil, is not significant for the much larger out-of-plane 
recoil. In addition, we note that, as seen in Fig. [6] the 
empirical formula underestimates the out-of-plane recoil 
to a moderate degree, which is consistent with the non- 
cancellation of the upward/downward recoils induced by 
the orbital plane precession. (We note that the amount 
of precession, per orbit, becomes progressively smaller for 
larger orbital separations, and hence the cumulative ef- 
fect of this non-cancellation is not expected to introduce 
a significant recoil at the orbital separations correspond- 
ing to the start of our simulations.) Thus it appears 
that there are higher-order effects in the precessing case 
that perturb the out-of-plane recoil as well. These newly 
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discovered effects are likely modeled by higher-order PN 
terms in the empirical formula and we are in the process 
of investigating how the in- and out-of-plane recoils are 
affected, with the goal of incorporating these effects into 
the formula. This will be the subject of a forthcoming 
paper by the authors. 

We also note that the that the difference in the magni- 
tude of the in-plane component of the recoil for the two 
choices of (r + , r , r_) are around 25%. These differences 
are due to the fact that the out-of-plane recoil is very 
large compared to the in-plane recoil, and hence a small 
error in the direction associated with the orbital plane 
will lead to relatively large changes in the calculated 'in- 
plane' recoil. For example, if the out-of-plane recoil is 
1000 kins" 1 and the direction of the orbital plane has an 
error of 5°, then the in-plane recoil will have an error of 
1000 sin(5°) kms" 1 ~ 90 krns" 1 . Note that this error is 
too small to explain the large magnitude of the in-plane 
recoil, which appears to be due to a newly discovered 
physical effect mentioned above. 



C. An alternative choice for the plane 

In this section we reanalyze the recoil velocities us- 
ing a choice of orbital plane adapted to each individual 
run. In Fig. [9] we show f(t) as a function of r during the 
merger. From the figure we can see that there is some 
scatter in the locations of the zero, maximum, and min- 
imum of r(t), but are generally independent of both the 
initial angle and mass ratio (this observation justifies our 
choice above of using fixed choices for (r + ,r ,r_)). In 
this analysis we chose r+ to be the location of the mini- 
mum in r(t), ro to be the location of the zero (i.e. location 
where \r(t)\ is maximized), and r_ to be the location of 
the maximum in r(t). This choice has the advantage 
that we do not use any 'fiducial' choices for (r + ,ro,r_). 
Table [V] summarizes the results for the transformed re- 
coil velocities and kicks, while Table [VTl summarizes the 
fit parameters A and B after fitting the results from 
each mass ratio to the form (jH)). A fit of the maxi- 
mum recoil velocity per mass ratio to the form §§§ yields 
H = 1.904 ± 0.041 when G = K and (H = 1.954 ± 0.05, 
G = (6151±804)kms- 1 ) when both H and G are allowed 
to vary. These results are in good agreement with those 
for the fixed choices of (r+,ro,r_). Finally, in Fig. [TD] 
we show fits of the maximum out-of-plane recoil versus 
mass ratio using this choice of adaptive plane parame- 
ters. Note that with this method the scatter in the plot 
is more significant, indicating that a fixed choice of plane 
parameters produces a more robust estimate of the out- 
of-plane recoil. We also note that the choice of plane 
parameters (r + , r , r ) = (2.2, 1.2, 0.7) gives the best fits, 
which is consistent with the observation that tq — 1.2 is 
well adapted to the locations of the maximum in |r| for 
the runs considered here. 

Finally, we note that the error in the recoil velocities 
due to finite starting time were of order 10 kms" 1 in 



TABLE V: The transformed spin (at r — ro), recoil veloci- 
ties in kms -1 , and angles between the in-plane QXXXTH000 
spins and the other QXXXTHYYY configurations in degrees 
when r+ is chosen to be the value of r when r(t) reaches its 
minimum, ro is the value of r when r(t) = (i.e. when \f(t)\ 
is at its maximum), and r_ is the value of r when r(t) reaches 
its maximum. Note that in the transformed system L z and S z 
are both negative (i.e. there is some partial spin/orbit align- 
ment. We denote quantities in the (transformed) orbital plane 
with a _L subscript. 

Config Vx Vji a±/m H a^/m" 



Q13TH000 52.1681 

Q13TH090 58.7557 

Q13TH130 213.996 

Q13TH210 37.8511 

Q13TH315 234.147 

Q17TH000 278.533 

Q17TH090 133.506 

Q17TH130 405.804 

Q17TH210 273.036 

Q17TH315 303.864 

Q25TH000 328.136 

Q25TH090 368.549 

Q25TH130 194.765 

Q25TH210 340.595 

Q25TH315 122.309 

Q33TH000 432.797 

Q33TH090 522.116 

Q33TH130 123.321 

Q33TH210 538.156 

Q33TH315 160.341 

Q40TH000 345.769 

Q40TH090 464.445 

Q40TH130 199.788 

Q40TH210 358.115 

Q40TH315 138.301 

Q50TH000 206.224 

Q50TH090 216.108 

Q50TH130 351.439 

Q50TH210 122.003 

Q50TH315 379.77 

Q66TH000 134.205 

Q66TH090 253.359 

Q66TH130 324.629 

Q66TH210 164.115 

Q66TH315 296.905 

Q100TH000 89.5687 

Q100TH090 119.496 

Q100TH130 52.4861 

Q100TH210 80.6623 

Q100TH315 60.4502 



-112.324 0.770573 
-101.697 0.7674 
-372.167 0.715948 
130.01 0.764949 
372.257 0.733002 

584.276 0.73713 
345.022 0.747727 
305.157 0.796759 
-577.037 0.729478 
-110.066 0.797954 

-965.098 0.74441 
-909.203 0.767864 
59.5598 0.777913 
951.677 0.758298 
-205.964 0.771452 

-1210.91 0.784273 
-865.374 0.797412 
392.695 0.764043 
1015.2 0.798003 
-603. 0.755614 

-1460.64 0.777363 
-1240.18 0.795298 
84.8549 0.779402 
1413.71 0.783135 
-323.804 0.77561 

30.7051 0.780222 
1263.93 0.77053 
1578.89 0.789517 
-618.748 0.768792 
-1493.21 0.797621 

890.868 0.769565 
1717.31 0.797964 
1042.76 0.800428 
-1364.25 0.769748 
-840.357 0.799556 

1423.6 0.793956 
1015.99 0.804421 
-244.07 0.789395 
-1485.43 0.803926 
387.967 0.788101 



-0.236828 0. 

-0.246703 1.29153 

-0.357787 -55.7042 

-0.244748 -178.98 

-0.319972 105.641 

-0.307434 0. 
-0.292225 56.9267 
-0.0790985 -66.0061 
-0.32118 -170.825 
-0.0710164 87.869 

-0.277851 0. 

-0.203131 -20.4251 

-0.186817 -111.668 

-0.24157 167.113 

-0.210448 73.1503 

-0.156319 0. 
-0.0611074 -36.3755 
-0.242325 -87.9496 
-0.0555484 148.237 
-0.263486 78.0326 

-0.199695 0. 
-0.0996155 -31.0583 
-0.184948 -97.0931 
-0.174788 177.471 
-0.207095 73.4494 

-0.185459 0. 
-0.230075 -47.8033 
-0.13106 -98.1525 
-0.235123 158.316 
-0.0947157 68.7318 

-0.233273 0. 
-0.106154 -50.8354 
-0.0583467 -100.477 
-0.233049 171.955 
-0.0682437 72.8232 

-0.16161 0. 
-0.0421833 -66.5996 

-0.14723 -119.948 
-0.0816034 151.372 
-0.157316 56.6673 
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FIG. 9: Md 2 r/dt 2 versus r/M as a function of mass ratio and 
initial spin orientation. Note that the locations of the zero, 
minima, and maxima are very similar. 
0.04 




-0.01 - 



TABLE VI: Fit parameters for wy = Acos((0 - B)tt/180) for 
each q. A is in units of kms -1 and B is in degrees for the 
adapted choice of plane parameters (r+,ro,r_). 



q 


A 


B 


1/8 


379.0 ±9.9 


108.51 ± 1.33 


1/6 


611.4 ±37.4 


353.47 ±3.33 


1/4 


948.5 ± 55.6 


169.16 ± 4.08 


1/3 


1353.4 ± 66.5 


195.32 ± 2.34 


1/2.5 


1427.5 ± 15.6 


176.70 ± 0.72 


1/2 


1620.5 ± 16.7 


271.46 ±0.60 


2/3 


1659.9 ± 82.8 


310.17 ±3.05 


1 


1504.9 ± 15.3 


340.60 ± 0.56 



FIG. 10: Fit of out-of-plane recoil (Vji/a) versus mass ratio q 
for the adaptive choice of (r+, ro, The top shaded region 
shows our predicted value for K = (6.0 ± 0.2) x 10 4 , the 
lower region shows the prediction based on the modification 
proposed by Baker et al., the lower (green) curve shows a fit 
to rf 1 , while the upper (blue) curve has a simultaneous fit 
to -q n and K. Note that the measured values of Vj|/a, while 
slightly overshooting the predictions, agree much better with 
our 0(r) 2 ) form. The dots are the locations of the data points. 

V/ffi 




each direction for all runs. We estimated these errors by 
looking for a systematic offset in the recoil velocity when 
plotted as a function of time (See Fig fTTj) . These errors, 
while not insignificant do not change the results of our 
analysis because they are the same order as the errors 
due to extrapolations to infinity and the errors in the fits 
of the recoil. 



r 



IV. CONCLUSION 

In this paper we explored the merger recoil from pro- 
cessing black-hole binaries with a larger spinning black 
hole (with initial spin a/m = 0.8) in quasi-circular orbit 
with a smaller non-spinning black hole. We introduced 
techniques to determine the normal to the orbital plane 
at merger and thus decompose the recoil into its in-plane 
and out-of plane components. However, there are impor- 
tant open questions about the accuracy of the determina- 
tion of the orbital plane, the spin direction, and spin mag- 
nitude. All these quantities are measured in the highly 
dynamical region around the merger. In this work we 



have introduced techniques to begin studying this prob- 
lem that will need to be refined. The issue of the in-plane 
recoil is particularly important, because large in-plane re- 
coils imply that our heuristic formula needs to be mod- 
ified for strongly precessing binaries. Nevertheless, our 
results, as seen in Fig. [8l indicate that the out-of-plane 
recoil has an 0(rj 2 ) rather than 0{rf) leading-order de- 
pendence on the symmetric mass ratio. This result, while 
agreeing with our prediction, appears to disagree with 
the recent work of Baker et al. [57]]. Thus additional 
work with new configurations may be needed in order 
to determine leading-order dependence, or indeed, if this 
dependence is a function of the configuration. 
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